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Abstract 



The statistical correlations between defects in the two dimensional 
complex Ginsburg-Landau model are studied in the defect-coarsening regime. 
In particular the defect- velocity probability distribution is determined and has 
the same high velocity tail found for the purely dissipative time-dependent 
Ginsburg-Landau (TDGL) model. The spiral arms of the defects lead to a 
very different behavior for the order parameter correlation function in the 
scaling regime compared to the results for the TDGL model. 
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Typeset using REVT^ 



I. INTRODUCTION 



We study here the statistical properties of a collection of point defects generated 
during the evolution of the two-dimensional complex Ginzburg-Landau equation (CGLE)&I. 
We will be interested in that portion of the parameter space where the COLE, driven by 
random initial conditions, has a regime of defect coarsening where the density of defects falls 
off with a power law in time. Our interest here is in the statistical properties of these defects 
and ultimately properties of the associated order parameter driven by the dynamics of the 
defects. Initially we will focus on the velocity distribution of the defects and the spatial 
correlations between defects. 

The approach developed here is based on the use of a set topological invariants which 
are apphcable to a large set of systems which generate defects as a part of an ordering process. 
In particular one is led to a clean expression for the velocity of the defect cores in terms of 
derivatives of the order parameter field evaluated at the core position. This approach not 
only allows one to investigate equations of motion obeyed by individual defects, but opens up 
the possibility of treating, the statistical properties of an ensemble of interacting defects. We 
have, from previous worlS in the area of phase-ordering kineticJl, the analog of the Maxwell 
velocity distribution for a collection of phase-ordering defects. 

In the defect coarsening regime for the CGLE, the defect density, n{t), scales as L^'^it) 
where L{t) is a characteristic length which grows with t, and t is the time of the evolution 
of the system starting with random initial conditions. In these circumstances, as shown in 
detail below, the defect velocity probability distribution is given, as in the purely dissipative 
time-dependent Ginzburg-Landau (TDGL) case, by 



m) = r(-) -I l + VVt;^ , (1) 




where the characteristic velocity v ~ is given explicitly below. Similarly the defect- 
defect equal-time correlation function has the same form (see below) as found for the TDGL 
case. In the case of the correlations between defect densities at different times, we find some 
rather weak deviations from the results in the TDGL case. 

These results for the statistical properties of the defects, inspire one, using ideas which 
have been successful for treating the TDGL case, to look at the order parameter correlations. 
In this case we find results quite different from the TDGL case. This is due both to the 
spiral arms and processional motion characteristic of defects in the CGLE. The spiral arms 
render order parameter correlations shorter in range, compared to the TDGL case, and the 
order parameter correlation function shows the behavior ^ TT?{t)W{r / L{t)). The precessional 
effects are predicted to be prominent in the two-time order parameter correlation function. 
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II. BACKGROUND 

The complex Ginzburg-Landau equation can be written§ in the form 

dt^ = bV^tlj + (l -u\tlj\^)^ (2) 



where ip is a. complex field, and b and u are complex parameters. For the appropriate set 
of parameters (choice of b and u) we find on quenching from an initially disordered state, 
that the CGLE generates a set of coarsening point defects. The characteristic distance 
between the defects increases with time due to the annihilation process between defects and 
antidefectJ^. 

For b and u real Eq.(|D reduces to the dissipative TDGL equation which is the most 
widely studied model for phase orderingi. If we set 6 = m = irj in Eq. and take rj large, 
we find, and after a simple gauge transformation, Eq.(|) reduces to 

- ^5*7/; = V'V^ + (1 - . (3) 

This equation, the nonlinear Schroedinger equation (NLSE)0, gives a highly idealized descrip- 
tion of the low temperature properties of a neutral superfiuid. Unlike the TDGL system, the 
NLSE supports several conserved quantities. In particular the quantity / d'^r |'?/'(r, t)p does 
not change with time. This model supports the same defects as the TDGL model, but the 
dynamics of the defects are quite different. Two oppositely charged vortices in the TDGL 
model move along the line connecting them toward annihilation. In the NLSE the same two 
vortices move at right angles to the line connecting them. 



III. DEFECT CHARGED DENSITY AND VELOCITY FIELDS 

The approach developed here allows for a direct connection between a set of field 
equations, like the complex Ginzburg-Landau equation satisfied by an order parameter field, 
and the equations of motion of the cores of a set of defects. It has only recently been 
understood, as discussed below, that these expressions for the defect velocity reduce to the 
same form as found in pattern forming studies using very different arguments. 

The approach developed here is motivated by addressing the question: What is the 
probability of finding a defect a distance r from an anti-defect? In work on phase ordering 
kinetics wei developed methods which are convenient for handling such questions. A mo- 
tivating factor was the realization that in treating statistical properties of defects one does 
not want to work with formal structures which require an explicit treatment in terms of the 
defect positions. This leads to problems of specification of initial conditions. Instead we 
looked for a way of implicitly finding the positions of the defects using the order parameter 
field ip itself. 

Let us consider the case of two dimensions where we have point defects. The case of 
line defects can also be treateclil0 using these ideas but will be discussed elsewhere. The 
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basic idea is that the positions of defects are located by the zerosliiryj^ of the order parameter 
field ip. Suppose, instead of the positions Yi{t) we want to write our description in terms 
of the zeros of ilj{r,t). It is not difficult to see that the defect charged density has the two 
representations 

N 

p(r, t) = 5(^(r, t))P(r, t) = J] g,5(r - r,(t)) (4) 

i=l 

where qi = V{Ti)/\V{ri)\ = ±1, and V^r) is the Jacobian associated with the change of 
variables from the set of defect positions to the field ip: 

V = ^e,,,,V,,i^V,,r (5) 

where we sum over the Hi, e^^^j 2-dimensional anti-symmetric tensor and summation 

over repeated indices is implied. For later reference, the unsigned defect density is given by 
= |p(r,t)|. 

For systems where only unit charges are present, p is the topological charge den- 
sity. Notice that qi is well defined even for system like classical fiuids where the circulation 
associated with a defect is not quantized. 

The dynamical implications of this approach are simple. If indeed topological charge 
is conserved then we would expect the charge density to obey a continuity equation. It was 
shown in Ref . ( ^) that p satisfies a continuity equation of the form 

dtp = -V ■ (pv) (6) 

where the defect velocity field v is given explicitly by 

^^a = -^Ee"/5(^V^^*-^*V;3^) . (7) 

where T> is defined by Eq.(^ and we must remember that v is multiplied by the defect core 
locating 5-function in p in Eq.(|]). Eq.(|^) gives one an explicit expression for the defect 
velocity field expressed in terms of derivatives of the order parameter. This expression for 
the defect velocity seems to be very general. Notice that we have not specified the form of 
the equation of motion for the order parameter only that the order parameter be complex 
and d = 2. For the CGLE our expression for the defect velocity reduces to 

= E (^v^^v^^* - b*v^r^p^) . (8) 

Does this expression for the velocity agree with our expectations for known cases? Let us 
assume that we have a defect of charge m at the origin of our two dimensional system and 
write the order parameter in the form: ip = i?e*^, R = r'^'e'" and 9 = mcj) + 6b^ where 
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again r and (j) are the cylindrical coordinates relative to the core at the origin. It is then a 
straightforward bit of calculus to show that the velocity given by Eq.® reduces to 



2b" I VJb + a E V/3u; ]-2b'iv^w-^J2 ^cpVpOs 1 • (9) 



\m\y / \ \m\ ^ 



If we ignore the contributions due to the variation in the amplitude, w, Eq.d) reduceli to 
Va = 2h"V oiOb + &'|^ Z]/3 ^ai3^ p^B- The first term is the only contribution in the NLSE case 
and states that a vortex moves with the local superfiuid velocityEj. The second term is the 
Peach- KoehlerS term first found in this context by Kawasaki0. These are the results from 
the phase-field approach and leads, for example, to the same type of interaction between two 
vortices as found in fluids. The velocity of a single isolated vortex is zero. For a set of two 
isolated vortices one has the expected behavior for the TDGL and NLSE cases. 

For our purposes here the more important point is to consider the work of Tornkvist 
and SchroderEl. Using methods of differential geometry, they looked at the derivation of the 
form of the velocity of a defect in the case of the CGLE. They comment, "The evolution of 
a system with (spiral) vortices may be described in terms of the defects, or filaments, along 
with values of the fields" R and 9 "at positions away from the defects of filaments. Such 
a separation into collective coordinates and field variables is non-trivial, and the present 
work comprises the first exact treatment of this kind for a dissipative system". The final 
equation they obtain, in our notation here and for two-dimensional systems, is precisely 
given by Eq.(||). Thus the velocity given by Eq.(|D reproduces the most sophisticated results 
obtained using other methods. 



IV. AUXILIARY FIELD METHOD 
A. Overview 

How can we use these expressions for p and ^{ip) to compute the measureable statisti- 
cal properties of an evolving CGLE system? We will use a generalization of an approximate 
method which has led to good results for the TDGL case. The basic idea is to assume that 
there is mapping from the order parameter field onto an auxiliary field m which shares the 
same zeros as if) in space. In particularly we require p[%l)\ = p[m] and = v[m] where 
again we use the result that the velocity is multiplied by a defect-zero- finding 5-function. 
These requirements are not very constraining since they only require that be proportional 
to m for small m with corrections that are cubic in m. It has been convenient to think of 
m(x) as a 2-vector whose magnitude gives the distance from x to a defect core. Thus, as 
discussed in more detail below, near the core we can take Re = rrix and Im ip = my. 

The main assumption! in the theory is that the field m is gaussian and the variance in 
m is determined by requiring that the defect charge density continuity equation be satisfied 
on average: 
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dti 
where 



dG,,{12) ^ ^ p(l)v(l)p(2) G,p(12) (10) 



G,,(12) =< p(l)p(2) > (11) 



and p(l) = p(ti,xi). With these assumptions and assumptions about the initial conditions, 
one can work out all of the statistical properties of the defects including G'pp( 12) and the 
defect velocity probability distribution function defined by 

nP{V) ^ {\p\5iV - vm , . (12) 

Corrections to this gaussian approximation can be investigated using methods of the type 



developed in Refs. 19,20 



The procedure then is to first compute Gpp{12) and Gjp( 12) assuming that m is 
a gaussian field. This will give Gpp{12) and Gjp{12) as functions of the auxiliary field 
correlation function 

C„;3(12) =< m,(l)m^(2) > (13) 

where a and /? take on the values x and y. Inserting these results for G'pp(12) and Gjp(12) 
back into Eq. (p!oD gives an equation for C(12). It will turn out that this equation for C(12) 
can be solved analytically. This result can then be fed back into the result for Gpp(12) to 
obtain an explicit expression for the defect density correlation function. As part of this 
calculation we obtain the average defect density n =< |p| >. Finally we can carry out the 
average over the gaussian variable m to obtain -P(V) as a function of C(12), and in turn 
obtain an explicit expression for -P(V). 



B. Expressing Gpp in terms of C 

The defect density in the defect-defect correlation function defined by Eq.(pJ]) can be 
written explicitly in terms of the gaussian auxiliary field m in the form 

P(l) = ^ep,^2e^,^2Vp,m^,(l)V^2m^2(l)5(m(l)) (14) 

and we sum over all the indices z/ and p. In the isotropic case, worked out previously, the 
evaluation of Gpp for the n-vector model for the general case of n = ci was facilitated by 
the decomposition of the average for Gpp into a product of averages corresponding to each 
component. This decomposition is not possible here because the complex coefficients in the 
CGLE couple the components of the order parameter as the system evolves. We need a more 
general approach. This more general approach involves using the general identity valid for 
Gaussian fields: 
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< m,(l)F[m] >= E / ^^2C?'^2 a.' (12) < — ^--F[m] > . (15) 

Using this result for all of the fields in Gpp acted upon by a gradient in Eq.(|l2), one can 
bring all of the gradients outside the average. This generates many terms which are products 
of the matrix C and averages proportional to the quantities 

g{l2) =< (5(m(l))(5(m(2)) > , (16) 

and 

and similar higher-order derivatives of the 5-functions which do not contribute to the final 
result. 

A key assumption in the evaluation of Gpp is that the system is isotropic in space and 
we can write: 

C,{12) = C.{rMh) , (18) 
V«a.'(12) = 0(12)r;. , (19) 

and 

V;.^V^^)a.'(12) = - [^,^.,(12) - Cj,,(12)] f,v - CJ.'(12)5,,' (20) 
where r = xi — X2, and 

C,V(12) = 0(12) (21) 

Cj,,(12) = iO(12) (22) 

and the primes in the superscripts indicate derivatives with respect to r. Using these results 
one can then carry our the sums over the spatial coordinate labels, the yu's in Gpp to obtain: 

G,,(12)=G«(12) + G(?(12) , (23) 

where 

GW(12) = ^(12)e.,.,e.;.,C,^^,,(12)Cj^,,(12) (24) 

and 
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It is easy to evaluate, using the results from Appendix A, the remaining averages over the 
auxiliary field: 

g,„,{i2) = r>"c„„,(i2)a(i2) (26) 

5(12) - ^ (27) 



where D is defined by Eq. (|A10 ). Expressing Cpi,i{12) in terms of Co(12) and A(12), as given 



by Eg . (|A2|) , and doing the sums over the i/'s, we find after some rearrangement, the result 
for the defect density correlation function: 

where 

7l = (l-/l)-^ (29) 



fT = V/o + (30) 



/o = , ^° (31) 



5o(l)5o(2) 



Ao = , ^ (32) 



where So{i) = Co{ii) and 



5o(l)5o(2) 



Q = {Q' + {A',f . (33) 



We still need to determine the auxiliary field correlation functions Co and A. It is easy to 

see that the result given by Eq.(^), in the isotropic limit where A = 0, reduces to result 
first reported by HalperinEJ 

^..(12) = II {h') (34) 

where 

h = ^ . (35) 

271 ^ ' 
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C. Satisfying Conservation of Topological Charge 



The calculation of the current contribution of Gjp on the right-hand side of Eq. (|T^) 
is much the same as for Gpp except for terms which involve the on-site correlation function 



S^'\l) = \< (Vm(l))^ >= - {y^CoirMti. 



r=0 



(36) 



Gjp is also proportional to the factor (2^rdr ^^^^ performing an integration over r, 

(37) 



we obtain the averaged conservation law, given by Eq.(|l^), can be rewritten as 

^ (Q7I) = 26'M + 2h"N 

where 

M = 7^Q (^o(l) + /oVVo + AoV^Ao) + 7I (/oVVo + A^V^A;, 

and 

N = i^Q (/oV^ Ao - AoVVo) + It (/o V'A[, + A[,VVo 
where we have introduced the time-dependent quantity 

^(2)(1 



^o(l) 



(vVo(r, 



r=0 



Eq.(p7D looks very complicated but simplifies if we replace /o and Aq with 

/o = fx cos n 

and 

Aq = /r sin fl . 

Then Eq.(P7|) can be rewritten as 



7t/t [It (2/t - R 
where 



+ {2 It -R)+ 1t^' [2/|f2' - f^S + hS' 



R = 26'cuo(l)/r + 2b' A + 2b" B 



S = -2b' B + 2b" A 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



and 
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(46) 



B = 2V/t ■ VQ + JtV^Q . (47) 
A solution to Eg. (^31) is given by 

2/t = R (48) 

and 

2flh' = f^S - hS' . (49) 
This last equation can be reduced to 

2fTh = -S (50) 



The set of coupled equation given by Eqs. (^81) and (|50D are equivalent to the equations for 
/o and Ao given by 

fo = b'{cuo{l) + V')fo + b"V'Ao (51) 

and 

Ao = 6'U(l) + V2) Ao-6"VVo . (52) 



This is the set of equations which must be solved self-consistently to obtain the unknown 
quantities /o and Aq and uj{l). 



D. Auxiliary Field Correlation Function 

Equations (j51|) and (|5^ ) are reduced to a set of differential equations in time if we 
Fourier transform in space and put in the time labels explicitly: 

d 

— /o(q, tit2) = a{q, ti)/o(q, ^1^2) - /?qAo(q, ^1^2) (53) 
^Ao(q, ht2) = a{q, ti)Ao(q, ^1^2) + /3Jo(q, ^1^2) (54) 

where 

a{q,h) = b'{uJo{ti)-q^) (55) 

and 
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/3, = 6V • (56) 
Equations ( ^3]) and (|5^ need to be solved together with the symmetry condition 

/(q, h, h) = /o(q, ti, t2) + ^ Ao(q, h, h) = /*(-q, ts, ^i) (57) 
and the initial condition 

f{q,to,to) = 2nfe-'^^'''^' ^g{q) . (58) 

This particular choice of initial conditions, corresponding to an initial correlation length i, 
is very convenient since all integrals can be carried out analytically for all times. Finally 
we must remember the normalization which follows from the definition of /(12) given by 
Eqs.(|3l|) and (H): 

/(II) = / (|^/(q'^i'^i) = l • (59) 
It is not difficult to construct the appropriate solution given by 

f{q,h,t2) = /?(ti,to)i?(t2,to)^7(g)e''''''(*^+*^-'*")e^^'(*-*^) (60) 

where 

R^t,,to) = e'^^o . (61) 

It is straightforward to take the inverse Fourier transform of Eq. (BOl) with the result 



/(12) = R{t^, to)R{h, to) e-K'^/^)' (62) 

where 

= f + Ah'T - 2ib"{ti - ta) (63) 

and 

T = ^i±^ • (64) 
We must stop here and satisfy the constraint given by Eq.(|59D. We have from Eq. (|6^) 

l = R'ituto)[Pj (65) 

where 
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L\h) = P{t,,h) = e^ + 4b'h . 

Eq.(|65D serves as an equation for co'o(ti) which can be easily solved to give 

2 2 

""'^^'^ = UK) ^ + ^b'h ■ 

Using Eq.(|65|) to express i?(ti,to) in terms of L(ti) we find 

/(12) = <|.(ti,t2)^^e-K'^/^r 
I — too 



where 



and 



L^(T) 



This last definition implies 



and 



e _iy2 



f{l2) = ^{hMh 

1 — lUJ 



where 



X 



and 



1 +CU2 

X = r/L{T) 



1 2 



There are a number of comments relevant to this result for /(12) given by Eq. 
First note that there is consistency between the definition of u!o{t) given by Eq.(^OD and 
solution for / which leads to Eq. (|67[) . For equal times, ti = t2 = t, we have 

/(r,t) = e-^-' 
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which is of the same form as in the purely dissipative caseEyi21l with a characteristic length 
L ^ VlTt. It we look at the on-site r = autocorrelation function, 



f{0,h,h) = ^tuh)^±^ , (77) 

1 + tu^ 



we can write for ti, t2 ^ ^o, 



Ao 



*(*!,*.) = (4^1 . (78) 



For ti ^ t2, uj approaches a constant and the nonequilibrium exponent Aq for $ also governs 
/(O, ^1,^2) and is given by Aq = 1 which is the same0 as for the TDGL case for n = d = 2. 

The main new result is that for non-equal times the auxiliary field correlation function 
shows a novel oscillatory behavior. One of our chief goals below is discuss the possibility 
of observing this phenomena. We note here that /(12) does obey a form of scaling for 
ti,t2 ^ No- 



where r = ti/t2, 



/(12) = fix, t) = $(r) / . ' e-^^TTZ^ (79) 

1 — tUJiT] 



$(r) = ^ (80) 



b" (T-\ 



and 

2 1 + LO-^Kj) 

Rather than discussing this result for /(12) in more detail, it is prudent to remember 
that /(12) is not itself directly observable. Thus let us turn back to observables and their 
dependence on /(12). We delay discussing the details of the oscillations in /(12) until after 
discussing how this feeds back into the determination of observables. 



V. DEFECT-DEFECT CORRELATION FUNCTION 
A. General Result 

Given the explicit solution for /(12), Eq . (|79|) , we can return to the evaluation of 
the density correlation function, G'pp(12), given by Eq.(^). The input we need for its 
determination is 7^^ = 1 — F^, where 
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and 



Q = ifof + i^of 



Inserting these results for 7^- and Q back into Eq.(^) gives 

<l>2(r) ( \ 

^'"'^^'^'^ = 2n^L^{T){l + u;^{Tf [{1 + u;^{r)) ) ^^^^ 

where 

e'(l - s) - 9'(t) 

In analyzing Gpp{12) we must be careful to distinguish the equal-time case from the unequal 
times case. 



B. The equal-time case 

If ti = ^2 = ^ and r = 1, the density correlation function can be written as: 

where 

, , e^'(l - 1 

9{x) = — — 



This is the same result found in the purely dissipative case. It is knownEl that the con- 
servation of topological charge for equal times requires one to include in the defect-defect 
correlation function the correlation of a defect with itself: 

Gppir,t) = 5ir)n{t) + Gppir,t) (89) 

where n{t) is the average defect density. Then conservation of topological charge is given by 

d^r Gpp{r,t) = . (90) 
14 



Inserting Eq.(p9D into Eg. (pOj) gives 

n{t) = - I d'rGpp^r^t) . (91) 



However using the form given by Eq.(p8l) we can do the integral in Eq. (|9l|) and obtain for 
the average defect density 

n{t) = \iui'^ = \— . (92) 

This is the expected result if scaling holds. One can also find n{t) by direct computation 
and obtain 

m = ^ (93) 



where uiq is defined by Eq.(^) and given in this approximation by Eq.(|67D. We see that the 
two determinations of n{t) agree. 



for Gpp(12) since 



C. The unequal-time case 

For the case r 7^ 1 we have that the conservation of topological charge holds directly 

/ dVG,,(12) = -lim^ = . (94) 

The final step follows since Q ^ for small r and 7|. is regular in this limit. If we set r = 
in Gpp( 12) given by Eq.(p5|) we obtain 



1 1 4r 

G'pp(0,ti,t2) = 2^,^,^y^ (1 + ^2(^))(1_^)2 ■ (95) 

We see that this quantity blows up a r — 1 signaling the existence of the 5-function at r = 
obtained for equal times. Thus we see that the limits r — and r — 1 do not commute. 



VI. DEFECT VELOCITY PROBABILITY DISTRIBUTION FUNCTION 

The defect velocity probability distribution function is defined by 

n(t)P(V) = (|p(^)|5(V - v[^])) = (|p(m)|5(V - v[m])) . (96) 

One of the main results from the last section is that at equal times the auxiliary field 
probability distribution is isotropic and has the same form as in the purely dissipative case. 
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This means that we obtain the same result here as found in Ref.( ^) and given earher by 
Eq.(|l]) where the characteristic velocity v{t) is given by 

v{t) = 2{b'f^ (97) 

where 

^W(t) = vV(r,t)|.=o-^o'W (98) 
and using the exphcit results for f{r,t), given by Eg . ([761) we obtain 



The result for P(V) given by Eq.(|I]) indicates that the probability of finding a defect 
with a large velocity decreases with time. However, since this distribution falls off only 
as V~'^ for large V only the first moment beyond the normalization integral exists. This 
seems to imply the existence of a source of large velocities. Assuming the large velocities of 
defects can be associated with the final collapse of a defect structure (defect-antidefect pair 
annihilation for point defects), BrayS used general scaling arguments to obtain the same 
large velocity tail given by Eq.(|l|). 

One probe of the defect dynamics is to study the correlation between two defects 
including the correlation between their velocities. The two-defect velocity probability distri- 
bution, P[Vi, V2,r], gives the probability that two defects separated by a distance r have 
velocities Vi, V2. This quantity was determined in Ref.( |23|) and since it is an equal time 
quantity the results found there hold here. The physical results from the calculationil of 
this quantity for the TDGL model, carried out in detail for n = d = 2 using the same 
approximations as indicated above, are relatively simple to state. For a given separation r , 
the most probable configuration corresponds, as expected, to a state with zero total velocity 
and a nonzero relative velocity only along the axis connecting the defects: Vi = — V2 = vx. 
Moreover there is a definite most probable nonzero value for v = Vmax for a given value of r. 
The most striking feature of these results is that for small r the most probable velocity goes 
as Vmax = 1^1 f and k = 2.19 in dimensionless units defined in Ref. 



VII. ORDER PARAMETER CORRELATION FUNCTION 

Thus far we have focussed on the statistical properties of the defects in the system 
and found results for Gpp and -P(V) which are very similar to the purely dissipative case. 
Only rather small differences arise when one looks at unequal times. Unfortunately neither 
of these quantities probes the full phase dependence of the auxiliary field correlation function 
which shows interesting oscillations in space at unequal times. We show here that this phase 



16 



dependence may be probed via the order parameter correlation function evaluated at unequal 
times. Indeed this quantity, within the approximate treatment given here, is quite different 
from the purely dissipative case even at equal times. 

The order parameter correlation function is defined by 

c^(i2) =< ri^mi) > (100) 

and our approach toward its evaluation will be to find the relationship between the order 
parameter and the auxiliary field m. In section IV we required that ip be proportional to m 
for small m near the core of a defect. In evaluating Eq. (|100|) we need a more general mapping. 
The procedure we will use here has been successful in the purely dissipative case@. Picking 
up on the point made in section IV, we choose m(x) to represent the distance from x to the 
closest defect. This physical picture can be realized by constructing ip{in) as a solution to 
the equation for a single stationary defect: 

bVl^ + {l + tuji - u\^\^) ^ = , (101) 

where m serves as the coordinate and ui is the oscillatory frequency. In the purely dissipative 
case, b real and ui = 0, one has for field points well away from a defect core: 

^(m)) = ^Poe'^^""^ (102) 

where, for a defect of charge n, 

</)(m) = n tan~^ {'my/rrix) . (103) 

In the purely dissipative case, insertion of Eq. ( |102|) for iplm) into Eq. (|10CI|) and carrying out 
the Gaussian average over m, leads to the resultE^: 

C412)=i,lfl^d.J1^0^ (104) 

with /(12) given by Eq.(B|) with a; = 0. This approximate result has been rather extensively 
tested in the TDGL caseE^ . 

In the CGLE case we have a new and interesting element. There is a range of 
parameters where one has spiral defects. Thus, unlike the TDGL case, one has spatial 
structure associated with individual defects beyond the core. In particular Hagan0 showed 
that the far-field solution of Eq. ( |101| ) is given by 



V^(m(x)) = ^oe*l'^W+9'"W) (105) 

where q is the wavenumber of the spiral arms asymptotically far from the defect core and 
V'o = i^oil)- Q depends on the particular parameters of the CGLE as discussed by Hagan. 
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While there are values for which q vanishes, as in the TDGL limit, we will assume that we 
work in a region of parameter space where q ^ 0. 

Using this set of mappings the order parameter correlation function is given by 

C412) = < e-*('^(2)+gm(2))gK<^(i)+«™(i)) > 
= f ^2^(1)^2^(2) e'i'l'^-^') e*9{-(i)-(2)) e-l^ (106) 



(2vr) 

where (Px{i) = x{i)dx{i)d(f){i) for i = 1,2. The action A and determinant D are given in the 
Appendix. In particular A is given by Eg. ( [A 1 3D in terms of polar coordinates and we have 
more explicitly, 

r)2 /.(X) POO 



0^12)=^',-^ x{l)dx{l) x{2)dx{2) e 

X e-5E,-'Ww^oWj(a;(i)^a;(2)) (107) 
where the angular integrations are given by 

J{x{l),x{2)) = r dm rd<f){2) e^Wi)-<^(2)) e«'^--(i)-(2) cosm)-m~o) (^Qg) 



and 9 is defined by tan 9 = A/Cq. Shifting the angular integrations we see that the 9- 
dependence factors out: 



J(x(l), x(2)) = 2tt e'^ r d(t) e'^ e^'^^"^^)"^^) cos0 _ (^gg) 

Jo 

If we change integration variables from x{i) to 

y^ = ^W^)x{i) , (110) 
we can rewrite Eq. (|107|) in the form 



C^(12) = ^o'7-7t' / yidyi / y2dy2 e-5(fi+^2)e^feyi-52y.) / ^frv.y^ cos^ (m) 

jT Jo Jo 2n Jo 

where 

q^ = q^/s^hT' ■ (112) 

Notice that the phase dependence of the auxiliary field correlation function is isolated in 
the overall factor of / in Eg. ( [L 1 1| ) . The integral over (p in Eg. (|1 1 1|) gives a modified Bessel 



function, but for our purposes we only need the power-series result: 
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1 r^TT ^ ^ 1 / fTV^V2\ 



2k+l 



2nJo " ^0 + 1)! V 2 / 



Etttt-^P^ (113) 



and 



2fc 

- A;!(A; + 1)! V 2 



C^(12) = ^o2(7^^f — 1— (^) J.(gi)J,*fe) (114) 



^fe(gi)= 1 ydye--ye''i^yy^''-^' . (115) 



where 



In the hmit g — » Eq. (|114|) does, after some manipulations, reduces to the result found in 
the TDGL case for n = 2. Note that for q ^ 0, except for r = |ri — very small, qi, 



given by Eq.( [112| ), is becoming increasingly large with yS'o(i) ~ L{ti). This means we need 
evaluate Jk{qi) only for large qi. Evaluation of Jk{qi) for large qi is facilitated by writing 

Jo 

where Jk{qi) = JkiQi, 1/2). We have then 

Jk{qi,a) = J{qi,a) (116) 

and 

J(g„a)= rdye-^y'e^'^^y 
Jo 



— e~3t +iJ"(g.,a) . 



(117) 



We see that the real part of J is exponentially small for large q'f. However, it is easy to see 
that for large qi 



, 1 2a 3(2a)2 
riqi, a) = - + — + + 



118) 



It Qi 



This means that the leading non-exponential contribution to the order-parameter correlation 
function comes from Jo(q'i) and is given to leading order by 



Jo(gi,a) 



^\ ,2a 
da) qf 



2i 



+ 



(119) 



Inserting this result back into Eq.( |114]) we obtain 
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C^(12)=V'o'(7^Vo(gi)Jo*te) 



12/ -2 



(120) 



2'^ (gig2)^ 

g6 (5o(l)^o(2))3/2 • 
The scaled portion of the OP correlation function for x 7^ can be written as 

W{x,T) = ^,sT{t{)Sf{t,)C^{l2) 
Wo 

= VlT ■ (121) 

For r 7^ 1, the correlation function depends strongly on rj = b" /b' via u. This is not true 
at equal times where 77 = b"/b' does not appear in the scaling function. For r 7^ 1 the 
oscillations in / are now clear in the order parameter correlation function. Writing out the 
real and imaginary parts we obtain 

2F 1 

W = — —2 [cos z + usin z] (122) 

V 1 + u;^ (1 — F^) 

2F 1 

W" = ^===- —^[-sinz + ujcosz] , (123) 

V 1 + uj'^ [1 — F^) 

where F is given by Eq.(|83D and u) and z by Eqs.(0) and (^2]). We are interested in the 
oscillations associated with 77 = b" /b' 7^ 0. These are most clearly manifested in W" and 
characterized by the zeros at 

uj = tan zq (124) 
The first zero as a function of scaled distance is given by 

xl = 2 + \u^ + ■■■ (125) 

for small uj and 

xl = 11^-2 + {l/uj) (126) 

for large u. 
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VIII. CONCLUSIONS 



By using some new ideas about how to characterize defect dynamics, we have shown 
how one can determine local expressions for the defect density and defect velocity in terms 
of derivatives of the order parameter fields. These exact results were then used to derive ap- 
proximate results for the defect-defect density correlation function, defect velocity probability 
distribution, and the order parameter correlation functions. Within these approximations, 
which work well for the purely dissipative case, we find that the results for the defect-defect 
density correlation function and the defect velocity probability distribution are substantially 
unchanged from the TDGL case. Thus these results seem robust. The results for the two- 
time auxiliary field correlation function indicate some interesting oscillation of its phase as a 
function of scaled distance. Since the defect-defect density correlation function depends only 
on the amplitude of the auxiliary field correlation function these oscillations are not present. 
In the last section we have seen that some remnant of these oscillations is present in the 
order parameter correlation function. However, another new element element for the order 
parameter correlation function is that the spiral arms for the defects render the interactions 
between different spatial points much shorter range than for the purel dissipative case. Thus 
for different spatial points at equal times the order parameter correlation function is down by 
a factor of n^{t) relative to the TDGL case. All of these results can be tested via numerical 
simulation. 

Acknowledgement: This work was supported in part by the MRSEC program of the 
National Science Foundation under Contract No. DMR-9808595. I thank Drs. I. Aranson 
and P. Kevrekidis for very useful discussions. 

APPENDIX A: 

In the purely dissipative case all correlations are isotropic: 

Cuu'iij) =< rn^{i)m^'{3) >= S^u'Co{ij) (Al) 

where = (1, 2). In the complex case, over time, the real and imaginary components of 
the order parameter are mixed, and this requires that we treat the more general correlation 
function for the auxiliary field 

Cuu'iij) = S,,u'Co{ij) + e^^/eijA(12) (A2) 

which satisfies the required symmetry for classical fields, 

a.'(u) = CMj^) (A3) 

if Co{ij) = Co{ji). Thus the variance of the gaussian field m is determined by the two 
independent functions Co(12) and A(12). 
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We will be concerned with various two-point averages over m of the general form: 
CABm =< A{m{l))B{m{2)) >= J d^x{l) d^x{2)A{x{l))B{^{2))^^{l), (x(2)) (A4) 
where the two-point probability distribution is given by 

$(x(l), (x(2)) =< 5(x(l) - m(l))5(x(2) - m(2)) > 



(27r)2 (27r)2 



w' ij 



1 1 

re 2^ 



{2nf{det C)V2' 
where the argument of the exponential is given by 



and the matrix W is the inverse of C defined by 



7fc 



W is given explicitly by 
where 



■ /5o(l)So(2) \ ■ 
"ij — zr7-\ ^ ^0 ~ ^0 



So{i) = Coin) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



D--" = So{l)So{2) - Cl 



(AlO) 



= Cl + 



and finally 



det C = D-^ . 



(All) 



(A12) 



If we express x(i) = Xi{cos 0i, sin the argument of the exponential in the distri- 
bution takes on the simple form 
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A^Y. ^?^o(0 - 2D-^CTXix-2Cos{(t)i -(f)2-0) (A13) 

i 

where 

Co = Ct cos e (A 14) 

A = Ct sin e (A15) 

and 
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